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q^ ; ABSTRACT 



We have studied the predictions for the LSP decay within the framework of a 

radiatively broken unified supergravity model without i?-parity. Assuming that 

qj ■ Higgs/slepton mixing is the only source of i?-parity breaking and responsible for 

[T i | the observed neutrino oscillations we obtain predictons for the LSP life-time and 

iy~. ■ branching fractions. 

1. Motivation 

| i S_e try B is Presently th e m0St P0PU1 , „ t „ S0lve th e „ 

problem of the standard model (SM). Here, the cancellation of quadratic divergences is 
guaranteed and, hence, any mass scale is stable under radiative corrections. The most 
economical candidate for a realistic model is the minimal supersymmetric extension of 
Q\ ■ the SM (MSSM). In the SM baryon and lepton number are protected by an accidental 

symmetry ( i.e. there is no gauge and Lorentz invariant term of dimension 4 or less 
that violates B or L via perturbative effects). This no longer holds in the MSSM 

P-i. due to the existence of superpartners. One way to assure Baryon and Lepton number 

conservation (and hence the stability of the proton) is to impose by hand a discrete, 

j> [ multiplicative symmetry called .R-parityH, R p = (— \y s + 3B + L ^ where S, B and L 

are the spin, baryon and lepton numbers, respectively. Aside from the long proton 
life-time, R p conserving models have the very attractive feature that the lightest 
supersymmetric particle (LSP) is stable and a good cold dark matter candidate^. 
On the other hand, there is strong experimental evidence for neutrino oscillations™™ 
which can be accounted for if R p is brokenQ 

In this paper, we will investigate the LSP life-time in a SUSY-GUT scenario where 
R p is broken explicitly via dimension 2 termsQo&t9. We have discussed this model 
detail in ref. [7] where the emphasis was on neutrino phenomenology in the frame- work 
of radiative electro-weak spontaneous symmetry breaking (REWSSB)lHl. Here, we 
are particularly interested in the implications for high energy collider phenomenology. 
We will focus attention of the case that the LSP is a neutralino. This is the most 
interesting case, since it occurs naturally over most of the SUSY parameter space. 
Note, however, that in models with broken R p there is no theoretical/cosmological 
prejudice concerning the color or electric charge of the LSP. The only requirement 
is that the LSP life-time is sufficiently short {tlsp ^ 1 sec) so as not to disturb big 



bang nucleosynthesis! 2 ! or sufficiently long [tlsp ^ 10 24 sec / B(LSP — > Xv e )] so as 
not to lead to an unacceptable distortion of the cosmic microwave background". 
The most general gauge invariant superpotential can be written as 

W = ^ Jk L I L J Et + y? jk L I Q j Dt-yf k HQ j U c k - f , I L I H+ (1) 

where the supermultiplets are denoted by a hat. The left-handed lepton supermulti- 
plets are denoted by Li (i = 1, 2, 3) and the Higgs superfield coupling to the down-type 
quarks is denoted by L . 

Let us first determine the meaning of the various terms of eq. |l|. Here, y§j k , y^ k 
and yV k denote the lepton, down-type and up-type Yukawa couplings, respectively, and 
yu-o is the Higgs mass parameter. However, in contrast to the SM the MSSM allows for 
renormalizable baryon [lepton] number violating interactions y® k [y^ k and /ij] . These 
couplings are constrained from above by experiment. The most model independent 
constraints can be obtained from collider experiment sEj or neutrino physics ta. It turns 
out that the individual lepton and baryon number violating couplings only have to be 
smaller than O(10 -3 ~ fewx 10 _1 ) with the exception of y® 2 \ ^ 10~ 7 from heavy nuclei 
decays. Thus, the R p violating couplings need not be much more suppressed than 
the lepton and baryon number preserving Yukawa couplings, (remember that e.g., 
y® lx ~ 3 x lCU 5 /cos/3). Somewhat stronger but more model dependent constraints 
can be derived from cosmologyE-3. 

However, the experimental exclusion area can be strongly enhanced by imposing 
theoretical constraints: in the minimal SU(5) SUSY-GUT model, the right-handed 
leptons, the right-handed up-type quarks and the left-handed quarks are embedded in 
a 10-dim representation, 10j = E^@Uf@Qi- The right-handed down-type quarks and 
the left-handed leptons are embedded in a 5-dim representation, 5, = D$ © Lj. The 
two Higgs doublets are embedded together with two proton decay mediating colored 
triplets, T and D , in 5-dim representations, 5o = -Do © L and 5 = T © H. Hence, 
both the lepton and the baryon number violating interactions arise from the term 

H-gut = -yijk^i'EjlOk , (2) 

where the boundary conditions at the GUT-scale, Mqut, are given by y^ k = yf k ^ = 
Vijk = Vijk- These relations, which predict the down-type quark masses correctly to 
within a factor of 3, are expected to also hold at a comparable level for the R p violating 
couplings. Thus, in general the baryon and lepton number violating couplings are 
correlated in SUSY GUT models. This leads to very strong constraints on any y^ k 
from proton stability which are much stronger than any constraint on individual 
Yukawa couplings^. However, it does not place any constraints on the coefficients of 
of dimension 2 terms, /i,, which are the subject of this paper. 

The outline of our paper is as follows: in section 2 we present the neutrino and 
sparticle spectrum obtained from REWSSB including /Zj. In section 3 we present 
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Fig. 1. Contours of constant mo, B$ and sin^ app /sin^5. in the tan/3-m!/2 plane. We set A$ = 
and [i = 2.5m 1 / 2 - 

the numerical analysis of the LSP life-time, tlsp, and LSP branching fractions. Our 
conclusions are presented in section 4. 



2. Radiative Electro- Weak Symmetry Breaking 

Without any assumptions based on theoretical prejudice there are many models 
with different SUSY particle spectra and vastly different phenomenology. Thus, it 
has become standard to derive the low energy particle spectrum from minimal su- 
pergravity model with only four independent parameters: the universal scalar mass 
parameter, m , the universal gaugino mass parameter, m^, and the universal A (B) 
parameter multiplying the tri-linear (bi-linear) terms in the superpotential [eq. IJ. 
This approach is supported by the observation that the absence of FCNC implies a 
high mass-degeneracy of all scalars with the same gauge quantum numbers (with the 
possible exceptions of the Higgs mass parameters). 

First, we have to minimize the Higgs potential given by 



V = (fx 2 + m 2 H )H^H+(^ixj + 
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where the low energy soft SUSY breaking parameters are obtained by renormalization 
group evolution below Mqut in the standard fashion. In order to stay as close to the 
notation of the MSSM as possible we follow our notation of ref. 
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v = y/vjv] , and tan (3 = v/v 



V2 ' l ~ V2 
and we parameterize the VEVs in terms of spherical coordinates 
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Fig. 2. Contours of constant (a) mj , (b) m~ b , (c) m^a (tree-level) and (d) m^o (incl. 1-loop radiative 
corrections) in the tan/3-mi/2 plane. The other SUGRA parameters are as in fig. |l|. 



Analogously, it is convenient to parameterize the R p breaking mass parameters in 
terms of three mixing angles 



tan # 3 = — , tan 6 2 = — , 

/i2 A* 1 cos (J 3 



tan 6*i = 



/A) 



COS #2 ' 



(6) 



and \x = JJJLiJJLi- The potential in eq. § is minimized by an iterative procedure using 
the analytic solution for tan#i = as our initial values. This procedure also works 
surprisingly well for tan6*i > 1. For small R p violating parameters we can also obtain 
very reliable analytic expressions in the basis where y\, is diagonal by expanding in 
powers of [li/fio 
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with the convention that indices in braces are not summed over. In general, one fixes 
the GUT input parameters are mo, vn\j2 and Aq. Sqo and /i are obtained by solving 
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Fig. 3. Contours of constant (a) m x o, (b) m ±, in the tan/?-mi/2 plane. The other SUGRA 
parameters are as in fig. fll. 

eq. [7| and [| while keeping tan j3 and v fixed. Here, we find it convenient to fix the 
fermionic spectrum given by \x rather than the scalar spectrum determined by rriQ. 

In fig. |l|(a) and (b) we present contours of constant GUT parameters m and 
B in the tan /3-mi/2 plane, and in fig. |](c) we show how well the approximation 
works for the minimization of the potential (eq. 0). We see that the deviation of 
the approximation obtained from eq. [^ denoted by sin 1 app from the true minimum 
obtained by numerical methods and denoted by sin 6[ is quite small as long as tan (3 = 
0(1 ~ 10). However, it breaks down for tan/3 > 40. 

2.1. Sparticle Spectrum 

From LEP experiments we know that there are no charged superpartners with 
mass below m z /2. Furthermore, we can deduce a similar constraint on the lightest 
neutralino which, in our model, is instable. In fig. |] and |3| we present contours of some 
relevant scalar and fermionic superpartner masses in the tan/3-mi/2 plane. We have 
chosen A = and \i = 2.5m 1 / 2 . The value of m is obtained by imposing REWSSB 
[see fig. 0(a)]. We see that the only relevant constraint arises from the experimental 
lower on the lightest neutralino mass denoted by M LSP . 



2.2. Neutrino Spectrum 

The LSP phenomenology of the model under investigation here is governed to a 
very good approximation by only one parameter, tan#i. This parameter also deter- 
mines the neutrino masses whose upper limits are given by 



m Ue < 4.35 eV , m v < 160 keV , m Vr < 23 MeV , 
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Contours of constant (a) m 2 , — m 2 , (b) sin 2 20 r „ in the tan /3-M^sp plane. We set A = 
= 2.5m]/2 and ?7ii/2 is replaced in favor of M^sp- 



While the is no direct evidence for non-zero neutrino masses, there is strong experi- 
mental evidence for neutrino oscillations which imply that the three neutrinos are non 
mass-degenerate. In this work, we will take the view that our R p violating terms are 
the only source of neutrino masses and, therefore, should account for all the existing 
neutrino mixing effects. Since not all experimental indication for neutrino oscillations 
appear compatible with each other or have the same statistical significance a selection 
has to be taken. a The solar neutrino puzzleP may be the most compelling evidence for 
neutrino mixing. However, the effect appears to involve only the first two neutrino 
flavors neither of which is likely to be the heaviest neutrino both on theoretical and 
experimental grounds. Thus, the solar neutrino puzzle is not well suited for a determi- 
nation of tan Q\. Instead we choose to solve the atmospheric neutrino problemB. This 



requires that we fix tan 6\ such that m„ 



m. 



10 eV and we set tan 8 2 = 1 (we 



use a small value for tan# 3 = 0.045 in order to solve the solar neutrino puzzle via the 
MSW-effectE3; this angle will turn out to be quite irrelevant otherwise). Over most of 
the parameter space under consideration here this implies m v <C m Vr ~ 0.1 eV. It is 
then straightforward to obtain lower limits on tlsp from upper limits on m„ T [eq. |l(J 
by simple scaling arguments. 



a It has recently been suggested that three neutrino flavor are enough to accomodate all three indi- 
cation for neutrino ocillationEIJ. However, for the sake of generality we will be more conservative. 
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In fig. ^ we present contours of constant values for m 2 u — ml and sin 2 26' Atl/T . We 
fix fi — ±2.5mi/2 and A = 0. We see that for positive values of \x the mass difference 
m l ~ r ViL, = O(10~ 9 eV 2 ) is very small and quite compatible with long wave oscillation 
(LWO)! 2 ^ solution to the solar neutrino problem [fig. f| does not change if we set 
tan# 3 = 0(1)]. For \x < there is a sizable region were m 2 — m 2 e = O(10 -5 eV 2 ) 
as required by the MSW explanation of the solar neutrino deficiencya. These results 
were first presented in ref. |^ 

3. LSP Phenomenology 

In this section, we will discuss the decay properties of the LSP. The main interest 
from the point of view of collider phenomenology is the question whether the LSP 
decays inside the detector (else the analysis is equivalent to the case of unbroken R p ). 

Since the magnitude of R p violation is parameterized by an priori free parameter 
tan#i, we cannot determine tlsp- The situation changes if we relate tan#i to the 
neutrino masses. In the first part of this section we will present the prediction for tlsp 
assuming the atmospheric neutrino puzzle is a result of R p violation. [This prediction 
can be turned into a lower limit by imposing any of the upper limits on m VT of eq. |T0| .] 
In the second part of this section, we will discuss the branching fractions of the LSP 
which is independent of tan#i and, hence, also of any assumption about neutrino 
masses. 

3.1. LSP life-time 

As pointed out in ref. 18] in models without R p the neutrinos and neutralinos are 



6 Note that fig. 7(a), fig. 8(a) and fig. 9 are mislabeled in ref. M. The region with \i > and the 
region with fi < should be interchanged. 
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indistinguishable. As a result, the formulae for the neutralino radiative decay in the 

ri in 

MSSM E.3 and the decay into three fermions e3 can be directly generalized to our 
model. However, we do have to include the effects of Yukawa couplings which were 
neglected for R p preserving three-body decays ^3. Our complete set of formulae will 
be presented elsewhere^. Here, we simply present our numerical results. 

In fig. [| we have plotted Tlsp vs. tan#i. We find that there are simple scaling 



relations if R p violation is sufficiently small (say tan#i ^ 0.1): 



Tlsp , m Vx oc tan 2 8 t 



[x 



e,n,r) 



(11) 



The parameters tan 62 and tan #3 which govern the neutrino oscillations have only a 
small impact on the LSP properties. In fig. |6| we present contours of tan#i (fixed 
by imposing m 2 ^ — m 2 T = 10~ 2 eV 2 ) and constant tlsp in the tan (3-Mlsp plane for 
fi = ±2.5mi/2- We find that that required range of the Rp violation is tan#i = 10~ x ' 
(x = 2 ~ 5) with the upper (lower) limit corresponding to small (large) values of 
tan/?. The corresponding range of the LSP life-time is ctlsp — 1 m ~ 0.1 mm 
(for Mlsp = 40 ~ 160 GeV) and can easily be determined at forthcoming collider 
experiments. Furthermore, we need larger R p violation (for fixed m Vr ) for \x < due 
to cancellation among tree-level and one-loop contributions. For tan (5 > 30 ~ 40 the 
one-loop contribution even dominates over the tree-level result. For M LS p ^ m z the 
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two-body decays dominate and differ from the full width only by a few % [fig. [7|]. 
(A partial analysis of the two-body decays has been performed previouslyEE) 

So far we have assumed exact universality at Mgut- However, it has been pointed 
out in ref. 2? that the evolution from the Planck scale M P to Mgut can already have 
a significant impact on the sparticle spectrum. This is particularly important in 
SO(10) based models were the Higgs and slepton universality is violated via gaugino 
effects, since the Higgs (slepton) fields belongs to a 10-dim (16-dim) representation 
(remember: below M G ut non-universal effects arise only from Yukawa couplings 
while the non- universal effects above M GUT arise from the gauge couplings). We 
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can accomodate this effect by modifying the boundary conditions at Mgut 
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where we typical expect Rh = 9aGC/T/(4vr) ln(MGUTA^p) — —0.1. 

In fig. H we see that the effect of non-universal terms is very significant (small) for 
small (large) values of tan/? were the down-typ Yukawa couplings are small (large). 

3.2. LSP Bracking Fraction 

So far we have used results from neutrino physics in order to eliminate the R p 
breaking parameters tan6*j {i = 1,2,3). However, to a good approximation this 
dependence drops out if we consider the branching fractions. In fig. |S] we present the 
branching fractions as a function of Mlsp for eight different channels. The dominant 
decay mode is into quarks [(a) and (b) are first two generations only; (e) is the third 
generation] with a strong enhancement into bb (e) for small tan/?. Invisible decay 
modes (c) are typically below 10% and the radiative decay (d) is insignificant. The 
leptonic decays into t + t~ (f), £ + £~ (£ = e,/i) (g) and i^r^ (h) is typically 0(10%). 
For Mlsp ~ 100 GeV the situation becomes much more transparent by considering 
the two-body decays [fig. |10J- Here, there are only three relevant channels with 
B(LSP -»■ W±tT) ~ 0.5 and B(LSP -► Z°u), B(LSP -»■ h°v) ~ 0.25. 
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4. Conclusions 

We have investigated the LSP phenomenology in supersymmetric models without 
lepton number conservation. In any model of this kind, lepton number is violated 
spontaneously via sneutrino VEVs as well as explicitly. Both effects are of the same 
order and have to be studied consistedly. Assuming that Higgs-sneutrino mixing is re- 
sponsible for the observed neutrino oscillations we find that the LSP decays inside the 
detector. The life-time can be determined over a large region of the SUSY parameter 
space. The branching fractions for all relevant decay modes are we presented. 
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